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The theoretical evaluation of major nuclear structure effects on the asymmetry of allowed Gamow- 
Teller /3-decay rates in light mirror nuclei is presented. The calculations are performed within the 
shell model, using empirical isospin-nonconserving interaction and realistic Woods-Saxon radial wave 
£C) ' functions. The revised treatment of p-shell nuclei is supplemented by systematic calculations for sd- 

shell nuclei and compared to experimental asymmetries when available. The results are important 
in connection with the possible existence of second-class currents in the weak interaction. 
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I. INTRODUCTION 

One of the possibilities to look for the existence of second-class currents in the weak interaction is through the 
study of differences in allowed Gamow- Teller (GT) /3-decay rates of mirror nuclei 0, 0, IS @ • 
Weak processes in nuclei at low energy are well described by the effective V — A interaction : 

ff V -^=%4^ + h.c, (1) 
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where jt (j M ) is the hadronic (leptonic) current, and Gf is the weak interaction coupling constant. The most general 
Lorentz covariant form for the hadronic current reads Jl=Vf l + A^ with 

Vf, = i$ p gvlt, + fljCV^V + igsk tJl tp n , (2) 



2M 

An = iip P SU7/i75 + 7^yo>75^ + igpk^lz 



ip n (3) 



standing for the vector and axial-vector parts. Here fc M is the momentum transferred, M is the nucleon mass, ij) v 
(if> n ) is the proton (neutron) field operator. The six terms are the main vector, weak magnetism, induced scalar in 
P)l and the main axial-vector, induced tensor and induced pseudoscalar in ©, with gy, gu, gs, gA, gr, gp being the 
corresponding coupling constants. 

The weak hadronic currents in Eqs. ©~© can be classified pj according to their transformation properties under 
the (/-parity operation, Q = C exp(«7r72), i.e. the product of charge conjugation (C) and rotation over 180° about the 
2-axis in isospin space. Under this transformation, the main vector term, weak magnetism in (|2J) and the main 
axial-vector and induced pseudoscalar terms in A^ are even {QV^Q^ = V£) and odd (QAj^Q^ = —A 1 ), respectively. 
They are usually called first-class currents. The induced scalar term in and induced tensor in A^ transform, under 
Q, with opposite sign with respect to first-class currents. They are usually referred to as second-class currents (sec). 
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The interest in the related weak form factors stem both from the desire of understanding them in the framework 
on quantum chromodynamics and in the search for the physics beyond the Standard Model [6| . 

Constraints on weak form factors come from two hypothesis that follow from approximate symmetries of the 
Standard Model : the conserved vector current (CVC) hypothesis and the partially conserved axial-vector current 
(PCAC) hypothesis. In the limit of exact isospin symmetry, the CVC hypothesis ensures that : i) the weak form 
factors gv{k 2 ) and guik 2 ) are related to the electromagnetic ones; ii) the induced scalar form factor gs(k) vanishes. 
Therefore, the latter term can contribute only at the level of isospin symmetry breaking. 

The PCAC hypothesis does not put any limits on the value of sec induced tensor form factor gT(k 2 ). Thus, the non- 
conservation of the Q symmetry (due for example to differences in u and d quark masses and to the electromagnetic 
interaction) will result in a non-zero contribution of gx- 

It is an experimental and theoretical challenge to place limits on the possible presence of the induced tensor term 
in the axial- vector current. The search for evidence of its existence is longstanding 0, H|- 

In particle physics, the evidence can be obtained, for example, from G-parity non-conserving decay branch, such as 
r — ► Tjirv J3 , or from the angular distribution analysis after the decay r — * lottv involving both first and second-class 
currents @|. 

The study of weak processes in nuclei offers various possibilities in this search. There are two limits on the value of 
gr obtained from muon capture experiments, which are consistent with the absence of second-class currents 0, ITo| . 

In (3 decay, one can determine possible contributions from the induced tensor term in two types of measurements 
0. The first one is through correlation experiments where one measures the correlation between nuclear spin and 
momentum of emitted /^-particle in the decay of oriented systems, or the correlation between the momentum of 
the emitted /3-particle and of the radiation from a daughter state. A recent re-measurement of /3-particles angular 
distribution in the decay of aligned 12 B and 12 N gives a non-zero value of gr ; 11]. These experiments require the 
knowledge of the weak magnetism coupling constant, which gives another uncertainty. Only in the last years an 
independent way to extract sec, based on the (3-a angular correlation in A = 8 system, has been proposed [h2. Il3[. 

Besides, the second-class currents can give rise to different rates of mirror /3±-decay transitions |]J . This phenomenon 
was actively explored in early seventies, stimulated by first experimental observations of large asymmetries in the ft- 
values in the decay of light nuclei [l4T ] . The experimental asymmetry deduced from measured halflives 

where (ft)± refer to the (3± decays in the mirror nuclei, is known at present to vary on average from 1% to several 
tens of percent, being about (5 ± 4)% on average Jl5]. From the theoretical point of view, the asymmetry 8 (0} can 
have two origins : i) the possible existence of the sec term in the axial-vector current, i.e. non-zero contribution of 
gr (S scc ); ii) the breaking of the only approximate isospin symmetry between the decaying nuclear states (S nucl ), i.e. 

S = S scc + S nucl . (5) 

Assuming the perfect mirror symmetry for the nuclear states and impulse approximation, the induced tensor current 
gives rise |2( to the asymmetry i5 scc proportional to the induced tensor constant and the sum of maximum energy 
releases in both mirror transitions, (Wq" + Wq). The more refined calculations 0, including the off mass-shell effects 
and meson exchange corrections which are crucial in many-body physics |l6j| , predict the asymmetry provided by the 
second-class currents to be 

6 scc = -4—J+^(\L-2()(W + + W -). (6) 
9A 3gA 

Here A and £ are two parameters related to the coupling constants for off-shell and mesonic corrections, and induced 
tensor coupling constants, J and L are ratios of the matrix elements of meson-exchange currents between nuclear 
states (see Ref. for more details). 

However, it has been shown [rH fisl [ToL l20l| . that the nuclear structure effects provide the principal contribution to 
the asymmetry, from 10% to 20% for different nuclei, making the extraction of the possible contribution from induced 
tensor term an extremely difficult task 0, 0|. 

All early estimations [13, Ha . fl9l l20l| of S nucL were carried out within the nuclear shell model for a few p-shell nuclei 
only (namely, for A = 8,9, 12, 13), since the proper treatment of the sd-shell nuclei was far outside available at that 
time computational power. The basic ingredients of the calculations were the Cohen-Kurath interaction |21| with 
empirical isospin-nonconserving (INC) corrections [2(| and the radial wave functions obtained with the Woods-Saxon 
(WS) potential. The results obtained have been used to extract limits on second-class currents 

To be precise, the main contribution to S n comes from the difference in the matrix elements of the GT operator 
(err) between nuclear states, due to isospin mixing and differences in the radial wave functions, as we extensively 
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discuss below. The calculations of resulting asymmetry in ft- values performed by different authors [TtI ITsl \l$L l2(i| 
vary significantly, either because of different choices in the parametrizations or because of the more and more refined 
shell model technique used by the subsequent authors. 

Apart from this, there are some contributions to S nucl due to the higher order effects, most of which are usually 
incorporated as corrections to the statistical rate function /, i.e. forbidden matrix elements, known induced currents, 
such as the induced pseudoscalar (see Ref. Q for a complete list). Contrary to the calculation of the main part of S, 
the evaluation of higher order corrections, using different nuclear wave functions |l7l 120]. gives quite similar values 
and in any case does not exceed a few percent. 

Since the mentioned above estimations of the GT transition rates in mirror nuclei, no other attempts to refine the 
calculations have been undertaken up to now, in spite of the essential improvement in the knowledge of the low-energy 
nuclear structure. Wider and more accurate experimental data on (3 decay [l5j are available nowadays, including 
sd-shell nuclei. In Fig. 1 we show the asymmetry extracted from the measured /3-decay lifetimes of the selected mirror 
transitions for most of which the asymmetry is different from zero by 2a or more. 

Besides, the theoretical description of p-shell and srf-shell nuclei has made much progress during the last thirty 
years j^. Since the limits on the induced tensor current are still not settled j^. [j| |^ ficl ITlj. it is worthwhile to 
improve the constraints on second-class currents coming from the /3 decay of nuclei using on one hand the current 
experimental systematics and on the other hand improved calculations. 

The aim of this paper is to present up-to-date shell model calculations, for p-shell and, for the first time, sd-shell 
nuclei, of the main nuclear structure effects on the asymmetry S nucl JSJ) arising from the difference of GT operator 
matrix elements due to isospin mixing in nuclear states and differences in the radial wave functions. The higher 
order corrections can present a subject for a future work. In this study, we use p-shell [2l| and sd-shell model wave 
functions |2^| obtained from the INC Hamiltonians (24[, in a sufficiently large model space easily accessible now, and 
realistic radial dependence, and we compare the results with the ex per imental data. 

The theoretical framework, we choose, follows the one of Towner [2(j and goes beyond it in several respects. First, 
we use the better fitted INC interaction to take into account isospin breaking, both without and within the parent 
state formalism, similarly to the calculations of Ormand and Brown |25j] for superallowed transitions. Second, we 
use a different parametrization of the WS potential. We carefully investigate the sensitivity of the asymmetry to 
the way the tail of the single-particle wave functions may be determined by fitting the single-particle energies to the 
experimental separation energies. In fact, there are two possible procedures : I) to adjust the energy of the last 
unoccupied single-particle state only, as in Ref. |19t 120] ; II) to fix all the energies of the single-particle states in the 
considered shell at the experimental separation energy, as it is done in Ref. 18] for A = 12 only. Strictly speaking, 
none of them is exact. However, we will argue that the latter is more consistent with the shell model, while the 
former can even lead to "unrealistic" results when the last unoccupied single-particle state is unbound. We will then 
use, contrary to previous calculations 0, |2(j , the procedure of fitting all the single-particle energies at the same 
experimental separation energy to estimate the asymmetry for an ensemble of p- and sd-shell nuclei. We fit the 
separation energy either by adding a surface term to the WS potential or by varying the volume term, as it is done 
in Ref. EHdi. 

In previous works on asymmetries and second-class currents [TtI ITsl IT^. I20] , it was never stressed that one should 
discuss the -B(GT)-values before considering the corresponding asymmetries. Here we carefully examine the B(GT)- 
values and test the wave functions by comparing certain electromagnetic observables, namely magnetic and quadrupole 
moments, to available experimental data. The problem of the effective GT transition operator in relation with the 
quenching of the GT strength is also discussed. 

We remark that to extract limits on second-class currents in (5 decay one needs calculations on B(GT)-values as 
precise as possible and therefore should concentrate on the best described nuclei. In the context of the shell model, it 
means that very light nuclei, possibly viewed as few-nucleon systems, may in principle not be very good candidates. 
For example, the lightest nuclei, such as those involved in the A = 8, 9 cases, should be treated with particular 
techniques. This is also why, contrary to what was done in the past, in the present work we make an extensive 
analysis of both p-shell and sci-shell nuclei (Fig. 1). 

In the present study we restrict ourselves to the WS wave functions. As was noticed in Ref. |25|, the WS potential 
may overestimate the difference in the tails of proton and neutron radial wave functions due to differences in the 
Coulomb and isovector potentials. The results shown in |25| indicate that calculations with the wave functions 
obtained from the Hartree-Fock procedure would lead to a smaller effect on the asymmetry. However, as we will show, 
the present level of accuracy on the J3(GT)-values is still not good enough for extracting possible contributions of 
the second-class currents from the asymmetries. In this respect, calculations with Hartree-Fock wave functions would 
not improve the -B(GT) estimates enough to change these conclusions. For these reasons, we present results obtained 
with WS wave functions only which provide the upper limit of theoretical asymmetries. 

The paper is organized as follows. In Section II we describe the theoretical approach, we give the details of the 
interaction used, the parameters of the single-particle WS potential and discuss the different procedures. Section 
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III is devoted to the presentation of the results on magnetic and quadrupole moments, transition probabilities and 
asymmetries obtained for p- and sd-shell nuclei. Summary and conclusions are given in the last section. 

II. THEORETICAL FRAMEWORK 

Within the nuclear shell model [2(|, the asymmetry Q is related to the GT matrix elements M± for (3 + and (3- 
decay through 



M_ 



M, 



(7) 



with 



M± = (f\\E n k=1 Hk)r±(kW) 



£ (-l)' / /+^+^+V(2J l + l)(2J/ + l) 



J f Ji 1 



ji h Jtt 

S 1 l\j2^J f J,)S 1 l 2 {j 1 -,W{n 2 l 2 j 2 \\aT ± \\n 1 l 1 j 1 ) , 



(8) 



where n is the number of valence particles in the many-body wave functions that describe the initial, and final, 
|/), states in mother and daughter nuclei with angular momenta Jj and Jf respectively; J^ is the angular momentum 
of the parent states n of the (A — 1) nucleus and S x / 2 {j; J J^) denotes a spectroscopic amplitude with a proper 
phase factor, and gives the mixing of different shell model configurations in the eigenfunction. The INC terms in the 
effective shell-model interaction are responsible for the difference in the values of S 1 ^ 2 ^; J J-*) for mirror transitions. 
The single-nuclcon reduced matrix clement is defined as 



(nihhWar+Wmhh) = <5 W2 (-l)^+ 3 / 2 V6V (2ji + l)(2j 2 + 1) <j 1/2 J2 h ) «• 

h 1/2 1 



JU2 



(9) 



with (m, k,ji) labelling the single-particle states in a spherical basis and ^J U2 being the overlaps of the single- nucleon 
radial wave functions : 



K^ h {r)Rl M2 {r)r 2 dT 



3132 I '"n 1 l 1 ji y J^-n 2 l2'32 
10 



(10) 



The radial single-particle wave functions Rn 1 i 1 j 1 ( r ) are obtained from a spherically symmetric single-particle potential. 
For realistic radial wave functions, obtained from WS or Hartree-Fock potential with Coulomb and other charge- 
dependent corrections, the integrals (|10() are non-identical for two mirror processes. The difference in the matrix 
elements of GT operator for (3 + and (3- -decays © is extremely sensitive to the asymptotics of the radial part of the 
single-nucleon wave functions which depends on the particle separation energy through 



R(r) oc exp 




2m\E-E„ 



(11) 



where E and E^ are the energies of the initial state of the A-nucleus and of the final state in the parent (A— l)-nucleus, 
respectively, and m is the reduced nucleon mass [Tsj . 

In case the number of the parent states to be taken into account in (A — l)-nucleus is too large, one is obliged to 
introduce different approximations. There exist two basic possibilities to deal with the problem: either to truncate 
the number of states, leaving only those which mostly contribute to the matrix element JHJ, and then scale the result, 
or to neglect completely the dependence on n and approximate all integrals by those corresponding to the (A — 1) 
nucleus being in its ground state. We have used the latter procedure. In this case, the expression JHJl reduces to 



M± = OBTD(ji, j 2 ] AJ = \){n 2 li3\\\oT ± \\mhh) 



(12) 



where OBTD(ji,_j2; AJ = 1) are one-body transition densities pBf. 
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We will present the results on -B(GT)-values obtained using both Eq. (|12fl and Eq. JSJ, where we have taken 
into account up to as many as 100 states of a parent nucleus. Provided this number is still too small for a good 
approximation (as it is for some sd-shell nuclei for which the dimensions of the eigenvalue problem reaches a few 
thousands), we have restricted ourselves only to the calculations according to Eq. I|12l) . 

To summarize, within the shell model the asymmetry in the decay of mirror nuclei can have two origins: i) 
inequivalent S 1 ^ 2 values in (JHJ) or OBTD's in Ijl2(l for /3+ and /3_ decays due to isospin breaking effects; ii) difference 
of the overlaps between the proton and neutron radial wave functions l|10|) for mirror transitions. (These two nuclear 
structure effects were referred to, in the past, as effects due to Coulomb and charge-dependent nuclear forces and 
binding energy phenomena, respectively 20].) The two factors can be estimated independently and therefore we will 
discuss them separately. 



A. Isospin corrections to configuration mixing 



The isospin symmetry in nuclei is only an approximate symmetry due to the presence of the Coulomb interac- 
tion and isospin non-conserving nuclear forces. Within the shell model, both effects can be well incorporated in a 
phenomenological way. In Ref. 01 1 the INC Hamiltonians have been derived by supplementing several empirical 
shell model interactions (e.g., Cohen-Kurath for p-shell nuclei, the interaction of Zuker and co-workers j2?J within 
the (pi/2 1/2) shell model space for nuclei around 16 O, the interaction of Wildenthal [2H for sd-shell nuclei) by 
extra terms including isospin non-conserving ones. These are isovector single-particle energies for a given shell, matrix 
elements of Coulomb force and of the isovector and isotensor parts of nucleon INC interaction, whose strengths were 
fitted to the parameters of the isobaric mass-multiplet equation. The interaction yields the shifts of the levels in 
mirror partners often within 100 keV precision compared to experimental ones. We adopt these INC interactions for 
all our calculations within the p, (^1/2^5/2*1/2) and sd configuration spaces. 



B. Radial wave functions 

We determine the single-particle radial wave function R n ij (r) of a nucleon participating in the (3 decay by solving 
the Schrodinger equation with a WS potential : 

V(r) = ~V ws f(r) - vJ-± i/(r)]l • s + V c h(r) (13) 
r ar 



where 



f(T) = — CTTT (14) 



h(r) 



[l + exp(^)] 

for r > Rq 



r 



with i?o = ro(A — l) 1 / 3 and V c = Ze 2 . The parameters are chosen following the Bohr-Mottelson parametrization 
[2|, namely r = 1.27 fm, a = 0.75 fm, V ws = V - V NZ * {N - Z)t z /A, V so = V u * V ws , V = 50.5 MeV, V NZ = 32 
MeV, Vis — 0.22. N, Z, A are the neutron, proton and mass numbers of the nucleus. We have checked that this 
parametrization reproduces well charge radii. The agreement for p-shell nuclei is within 5%, except for the halo 
nucleus 8 B. The WS potential used in Ref. [3 does not contain any symmetry term, while the one of Ref. |20| differs 
in the parametrization and includes an additional Coulomb term. 

Remark that the potentials we use for the mirror partners differ only for the Coulomb term. 

In general, there are two possible procedures to determine the radial part of the wave functions by adjusting the 
energies of the single-particle states in the potential of the parent (A — l)-nucleus to the experimental proton or 
neutron separation energy in the A-nucleus. Method I : we fit the last unoccupied single-particle state and keep the 
same potential to get all the other radial wave functions. This choice was followed in Refs. fl9t I20I . Method II : we 
adjust each single-particle energy at the separation energy of the nucleus, as was done in Ref. jT^] for the A — 12 case 
only. We believe that the latter procedure is more consistent with the shell model than the former because method 
I corresponds to the extreme independent particle model. This conclusion is also supported by the unrealistic values 
for the electromagnetic moments in A = 8, 9 obtained with method I. 
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In practice, the adjustment of the energies is performed either by including a surface term (~ [f'(r)] 2 ) or by varying 
the depth of the volume term in the potential. We use the same set of parameters fitted for a given nucleus for its 
mirror partner. This automatically yields single-particle energies very close to the experimental separation energies 
(the difference does not exceed a few hundred keV and in some cases 1 MeV). 

We have studied the sensitivity of the asymmetries to both methods I and II (either including the surface term or 
by varying the depth of the volume term) for p-shell nuclei as trial case, and we have applied method II only, the 
most consistent one, to calculate the asymmetries for the se?-shell nuclei. 

III. RESULTS 

In this section we present the calculated magnetic and quadrupole moments, the transition rates and the asymmetries 
for an ensemble of p-shell and sd-shell nuclei. The S 1 ^ 2 values and OBTD's have been obtained from the shell model 
code Oxbash |2^. We use different types of radial wave functions : Harmonic Oscillator (HO) and Woods-Saxon 
(WS) wave functions obtained from the different fitting procedures as explained in the previous section (method I and 
II), and either neglecting the dependence on the parent state <|12[) or including it l|8 HlO[l . We also discuss how the 
results vary when one adjusts the single-particle energies either by including a surface (S) or by changing the volume 
(V) term of the WS potential. 

The experimental J5(GT)-values are obtained following the procedure outlined in Ref. [sol l8ll| with updated data 
on Q-values, half-lives and branching ratios as specified in Tables II and III (from the compilation of asymmetries of 
allowed transitions (log(/i) < 6, 7Tj = 7r/) in A < 40 nuclei presented in Ref. ^^|). For the decay to unbound states 
in the A=8,9,17 systems, we use the approximation of widthless levels and therefore we do not adopt the f- factors of 
Ref. HU. 

For comparison with theory, we have selected only the transitions for which the error on the asymmetry is small 
enough, since our final purpose is to extract the contribution of the second-class currents, if any. 

Before discussing the results on B(GT)-values, we comment on the problem of quenching of the GT strength. We 
also show what level of accuracy of the nuclear wave functions we have, by considering the electromagnetic properties 
of some of the nuclear states involved in the studied /3-decay branches. 

A. Effective GT-operator 

It is was suggested long ago that the axial- vector coupling constant qA in GT transitions is quenched in the 
nuclear medium compared to the free space value (for review see, e.g., Ref. j40| and references therein). In fact, the 
experimentally measured total GT strength from charge-exchange reactions is largely suppressed compared to the 
model-independent sum rule [4pl | , as well as the /3-decay transition rates are systematically overestimated by the Ohu 
shell-model calculations (e.g., Ref. 30] ). The reasons for this discrepancy are (i) lack of high-order configuration 
mixing, and (ii) inadequacy of the one-body GT operator and need of taking into account meson-exchange currents, 
in particular, A-isobar admixtures (the former effect is estimated to provide a dominant contribution to the quenching 
within the Ohui shell model |4lLl4^h 

However, it is also well established [3(], EH that in the shell model, up to a good approximation, the phenomenon 
can be taken into account by an empirical multiplicative factor, (gA/gv) e g — A {9A/gv)\ >alc with a bare value for the 
ratio of 1.26 as obtained from neutron decay (see Ref. 01 f° r certain cautions). In this case, the quenching does not 
influence asymmetries of the transition rates, i.e. the ratios of the matrix elements, which is then simply given by Q. 

Note, that here we are concerned with individual transitions, which represent only a small fraction of the total GT 
strength, even being the main branches of the p decay. Their analysis is more challenging because of the extreme 
sensitivity to configuration mixing predictions, in particular, for weak transitions or in the presence of close-lying final 
states. Thus, the quenching phenomenon as described above can be more clearly seen for stronger transitions. 

Bearing these remarks in mind, in the sections below we give the theoretical values for the GT-transition rates 
obtained with (gA/gv) halc = 1-26. 

B. Test of the wave functions 

The accuracy of the wave functions is one of the most critical issues of the study, especially since we are interested 
in the fine details of the weak interaction. To check the wave functions of the states involved in the (3 decay considered 
below and to demonstrate the present level of accuracy, we have computed the transition probabilities of electromag- 
netic processes for some of the states of interest. The electromagnetic processes are described up to an uncertainty 
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arising mainly from two factors, provided the best wave functions of the many-body problem are available: i) the 
restricted model space used in the calculations; and ii) impulse approximations for the operators. 

We have calculated some electromagnetic characteristics, namely, magnetic dipole and electric quadrupole moments 
for the ground states, which provide the most stringent test for the wave function of a given state. We will follow 
the conventional way and use the standard form of the Ml operator with the free-nucleon ^-factors (<7 s (7r) — 5.586, 
<?;(7r) = 1.0, g s {v) — —3.826, gi{v) — 0.0), which gives in general a reasonable agreement with the data. The electric 
quadrupole moments have been calculated with three different choices of effective charges [22, • The results are 
summarized in Table I and Fig. 2(a,b), where the calculated values for certain states in p-shell nuclei exploited below 
in the /?-decay analysis are compared with the data. 

Both types of moments were computed according to 

M(L) = OBTT)(j 1 ,j 2 ;AJ = L)(n 2 l 2 j 2 \\d L \\n 1 hj 1 ) , (15) 

31 ,32 

where Ol is magnetic dipole or electric quadrupole operator. The OBTD values have been obtained from INC 
interaction, while the radial wave functions from either the harmonic oscillator or the WS potential, adjusted to give 
the separation energies only with respect to the ground state. 

The magnetic dipole moments are sensitive to the single-particle configurations and almost do not depend on the 
radial form of the basis wave functions. As seen from Table I and Fig. 2, the theoretical description of the magnetic 
moments is always within 20%, thus showing the high quality of the effective interaction we use here. 

The quadrupole moments depend on both the configurations involved and the radial wave functions. The results 
obtained with both harmonic oscillator and WS radial wave functions (method II, V) for three choices of the effective 
electric charges are given in Table I. The calculations obtained by fitting the surface term (S) gives very similar results 
and to avoid repetition we do not present them. 

To get the radial wave functions, we have also tested the other possibility of adjusting the potential, namely, fitting 
the separation energies to the last unoccupied nucleon state (method I, V or S). We have found out that this procedure 
yields unrealistic values of the quadrupole moments for light nuclei. For example, the moments of 8 Li calculated with 
the radial wave functions I V and I S are —0.85 e.fm 2 and —0.88 e.fm 2 respectively, to be compared with the 
experimental values and theoretical results obtained with method II (Table I). We will see that the calculation of 
/3-decay transition rates does not help in choosing between method I and II, because they both yield similar agreement 
with the measured £?(GT)-values. We emphasize that only method I (V or S) was used in the past 0,|2(j, the case 
A = 12 being the only exception [l8j. 

As seen from Table I, both harmonic oscillator and WS radial wave functions (method II) yield similar values of 
the quadrupole moments and on average reproduce fairly well the experimental data. It is difficult to use the results 
as an additional constraint for the potential, since the quadrupole moments are also rather sensitive to the choice of 
the effective charges (see Table I) . The values obtained with the WS wave functions and the last set of the effective 
charges (e(7r) — 1.35, e{y) — 0.35) are plotted in Fig. 2(b) and as seen agree within 20% with the results of the 
measurements. 



C. GT-transition rates 

The calculated S(GT)-values, 

B(GT±) = Wrr |M±|2 ' (16) 

are summarized in Tables II, for p-shell nuclei, and III, for sc?-shell nuclei, in comparison with the experimental 
data (Fig. 3). The theoretical results are organized into columns labelled by the type of calculations. Contrary to 
the moments discussed above, the /3-decay rates are extremely sensitive to both, the asymptotics of the radial wave 
functions and to the correct accounting for the dependence on the parent state. First three columns with the results 
on B(GT) values in both tables are obtained using formula (|12|l . while in the last column (denoted by a star) the 
radial dependence on the parent state according to (jHJ) is taken into account. The label IC(INC) refers to the isospin- 
conserving (isospin-nonconserving) interaction used, and HO or WS indicate the type of the wave functions. Roman 
numbers I or II in Table II stand for the method to fit the single-particle potentials as discussed in section II. Only 
method II is used in the calculations for sd-shell nuclei. 
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1. p-shell nuclei 

A=8,9. 

The analysis of light nuclei is often complicated from both theoretical and experimental sides. The studied A = 8 
triad includes the transition to the very broad 2 + state of 8 Be, unstable with respect to act breakup. Moreover, in 
the description of the well-known proton-halo nucleus 8 B, the vicinity of the continuum can be important. The decay 
of proton-rich 9 C always continues through 2a-p three-body break-up, which complicates the experimental analysis 
of corresponding transition s l34l . In addition, the specific feature of 9 Li- 9 C mirror pair is the so-called anomalous 
isoscalar magnetic moment |48| which is still not fully understood from a theoretical point of view. 

The shell model yields good agreement of the energy spectra and even nuclear moments with either harmonic 
oscillator or WS wave functions (see Table I) for the lowest states of A = 8,9, despite missing the width of the 
unbound states. The only significant deviation is the quadrupole moment of the ground state of 8 B, showing necessity 
to go beyond the standard shell model we use. 

However, as is clearly seen from Table II, the shell model experiences difficulties in reproducing the absolute values 
of the GT /3-decay rates which are about ten and five times overestimated for A = 8 and A = 9 systems, respectively. 
The calculations with the best effective GT operator of Ref. [3lJ also overestimates the transition rates in A = 8,9, 
although with a factor of two discrepancy. 

The apparent disagreement is unlikely to be the consequence of a particular renormalization of the GT operator for 
these nuclei, and is an indication of the limits of the theoretical approach used (few-nucleon systems, clusterization 
effects, close continuum). Following these arguments, one should not consider data on A = 8,9 to extract second- 
class currents, until an accurate treatment by more fundamental microscopic techniques, such as Green's Function 
Monte-Carlo methods |52(, and no-core shell model [53j is attained. Note, that A = 8 nuclei treated by the standard 
shell model have mainly been used in the past in the context of the second-class currents problem (e.g., 

A=12. 

The triad of A = 12 nuclei, 12 B(1+;1) -> 12 C(0+;0) and 12 N(1+;1) -> 12 C(0+;0), is the most famous and well 
studied case |54| , since it was the first measured case of large asymmetry and up to now with the best precision. The 
electromagnetic moments for 1 + ,T = 1 states in mother nuclei are well reproduced. The rates of the GT transitions 
to the ground state of 12 C are close to the measured ones, while the calculated /3-decay rates to the first excited state 
are about two times larger than the experimental values. The use of effective GT-operator [3l| helps to reduce the 
disagreement to only 20% of overestimation. 

At present, there are already afr-initio shell-model calculations available for 12 C 55]. It would be interesting to 
apply the model to calculate the GT strength, since the questions such as limited configuration space and necessity 
of auxiliary realistic potential are not encountered there. 

A=13. 

The data on the rates for the GT transitions to the ground states (^ , 4) in A = 13 nuclei is also advantageous by 
its precision. The shell model well reproduces the four lowest excited states of negative parity, including the values 
of the ground state magnetic and quadrupole moments. The positive parity states are outside of the p-shell model 
space. The theoretical /3-decay rates are a factor of 1.2 larger than the experimental ones. No particular differences 
between wave functions obtained with the fitting procedure I and II can be seen in this case. 

2. sd-shell nuclei 

A=17, 

We will devote a special discussion to A = 17 nuclei, since they have already got a lot of attention from both exper- 
imental and theoretical sides, because of the enormously large asymmetries of the first forbidden /3-decay transition 
rates pffi 1571 158| . 

Here we study these nuclei with respect to the asymmetries in the allowed GT decay modes, which are also large. 
Although the nuclei near the doubly magic ones are the most difficult for the shell model framework, nuclei around 
16 can be well described within pif 2 d^/2 s i/2~^ ae ^ model space |27|. In the calculations we used the interaction of 
Ref. [24| , which is an INC version of the interaction close to the pioneer one |27j ■ Let us remark that the interaction 
exploited here gives much more realistic spectra, in particular the positions of 3/2~ states in 17 and 17 F, and 
transition rates as compared to the (psd)-shell model interaction used in Ref. [3l| . 

A=20,25,28,31,35. 

As is well-known, the sd-shell nuclei are fairly good described by the interaction of Wildenthal [2^, , while the 
INC corrections take care about the subtle effects. The systematic calculations of /3-decay rates in sd-shell nuclei 
have been already performed [22l l30|. although without accounting for INC effects, and are known to be in a good 
agreement with the experiment. 



9 



The results on _B(GT)-values for sd-shell nuclei, obtained with both IC and INC Hamiltonians and with harmonic 
oscillator or WS radial wave functions are summarized in Table III. There are still no data on the transitions for 
A = 21 and A = 25 nuclei, however, the measurements are under way j^. The calculated B(GT) values are all about 
a factor of 0.5 — 0.7 quenched compared to the experiment, except for the A = 20 case (T = 2 — > T = 1 transitions). 

As in the case of p-shcll nuclei, there was no principal difference between the radial wave functions, however, 
the values obtained with formula JBJ are assumed to be the most realistic ones. We did not take into account the 
dependence on the parent states for A = 25, 28, 31 because of the extremely large dimensions of the model spaces, as 
well as for A = 35 (S-CL-K-Ar) for which the experimental error bars are very large. 

D. Asymmetries of GT transition rates 

We have performed, for the first time, a detailed study of the asymmetries for a selected ensemble of p- and sd- 
shell nuclei. The chosen nuclei are on one hand rather well described by the shell model and on the other hand the 
corresponding measured asymmetries present small error bars (Table IV and Fig. 4). 

We have calculated the asymmetries of the B(GT)-values in the mirror nuclei due either to the use of the isospin 
non-conserving (INC) interaction or/and to differences in the tails of wave functions of the decaying proton or neutron. 
We present the effect of the isospin non-conserving interaction alone (INC+HO) , of the mismatch of the radial integrals 
alone (IC+WS) and of both (INC+WS). In the latter case, we show the influence of the inclusion of the dependence 
on the parent states in the radial wave functions ((INC+WS)*). 

We also present the first estimates of the asymmetries for some sd-shell nuclei for which the shell model description 
is rather good, but the asymmetries either present too large error bars or their analysis is under study [15| (Table V 
and Fig. 4). 

All the results shown in Tables IV and V correspond to the radial wave functions obtained with method II by 
varying the volume term of the potential in the fitting procedure. We have checked that fitting the single-particle 
energies with the inclusion of a surface term changes the results by a few percent. We have also calculated the 
asymmetries using method I. This always leads to much larger asymmetries than those determined by method II and, 
in particular, it gives unrealistic ones for cases where one single-particle state is unbound, such as for A = 8, 9. We 
have also mentioned above that it leads to unrealistic results for the quadrupole moments of some light nuclei. For 
these reasons and because of a better consistency with the shell model, we believe that method II (and not method 
I) should be used in the procedure of fitting the wave functions contrary to the previous studies [l^ . l2(i| . 

As we can see from Table IV, the effect due to the differences in the tails of the wave functions is by far the largest 
in the p-shell nuclei; whereas in the case of the sii-shell nuclei, it can be sometimes of the same order as the asymmetry 
due to the INC interaction. This happens when there is a large contribution of the d 5 / 2 or d 3 / 2 single-particle states 
for which the centrifugal barrier tends to confine the wave function inside the potential reducing the differences in 
the tails. The case A = 17, including the halo nucleus 17 F, represents an exception : the asymmetry due to the INC 
interaction is almost ten times larger than the one due to the mismatch of the wave functions. 

Concerning the dependence on the parent states (Eq. JSJ), we see from Table IV that its inclusion has a significant 
effect on the asymmetries. As it was found in 0,|2(]], the use of Eq.© instead of Eq. l|10[) may significantly increase 
or decrease the asymmetry and in a few cases up to a factor of two. 

We have compared the theoretical asymmetries to the ones obtained in Ref . 0, |2(J , although those calculations 
differ from ours in the parametrization of the shell model Hamiltonian, the WS potential, as well as in the fit ting 
procedure (method I was used). In general, our asymmetries are much larger than previously calculated in Ref. |l9Ll2Cl . 
However, for A = 12 case, our result is in reasonable agreement with their values, as well as with the value in |l8j 
where method II was used and also a very similar parametrization of the WS potential. 

As far as the comparison with the measured asymmetries is concerned (Table IV), we see that our calculated 
asymmetries due to the only approximate isospin symmetry in the nucleus (described in our case through the use of 
the INC interaction [24| and of WS radial wave functions) : a) are compatible with the measured asymmetries in a 
few cases; b) overestimate the measured values for A = 8, 12(g.s), 17; c) underestimate the experiments for A = 9, 20 
(Fig. 4). Remark that A — 8,17 represent rather delicate cases because they include halo nuclei which should be 
treated with particular techniques. For A = 8, 9, the theoretical £?(GT)-values are larger than the experimental ones 
by a factor of four. 

For p-shell nuclei, these conclusions do not change even after the inclusion of the other above mentioned higher 
order corrections to S nucl , which have been found in Ref. [2(| not to exceed 1%. 
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IV. CONCLUSIONS 



Within the shell model, we have calculated the Gamow- Teller f3 + an d P- decay rates for an ensemble of light p- and 
sd-shell mirror nuclei. The theoretical results are important in connection with the possible existence of second-class 
currents in the weak interaction. To assess their existence, it is necessary to carry out accurate estimates of the 
asymmetries coming from nuclear structure aspects. 

The calculations are performed systematically for the known experimental cases among p-shell and, for the first 
time, sd-shell nuclei. We have evaluated the asymmetry due to the only approximate isospin symmetry of the shell 
model Hamiltonian and differences between the proton and neutron radial wave functions. The former is accounted for 
through the use of an isospin non-conserving interaction which describes well the shifts in energy levels of the mirror 
nuclei. The latter is due to the use of realistic single-particle radial wave functions from a WS potential adjusted to 
reproduce the experimental proton and neutron separation energies. 

We have thoroughly studied the sensitivity of the asymmetries to the inclusion of the parent state dependence as 
well as to different fitting procedures, namely by fitting cither the energy of the last unoccupied state or the energies 
of all the single-particle states in the shell at the same experimental separation energy. We argue that, contrary to 
what was done in the past, the latter fitting procedure is more consistent than the former. 

The theoretical asymmetries we obtain are on average larger than all previous estimates. The results are compatible 
with the measured values in a few cases. For A = 8, 12, 17, the theoretical values are larger than the experimental 
ones, whereas for A — 9, 20 we find the opposite. 

As is clear from the present study, the average reliability of £?(GT)-values, which determines the level of confidence 
for the theoretical asymmetries, seems still not to be high enough for extracting the second-class currents from the 
differences in nuclear GT /3-decay rates, even at the present level of many-body techniques. 

In this work, we restricted ourselves to WS wave functions. Calculations with Hartree-Fock wave functions would 
not modify the quality of the B(GT)-values significantly enough to change these conclusions. Nevertheless, these 
calculations would be needed, once a better overall description of the transition rates is attained. 

More extensive and systematic data would be very helpful. 
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TABLE I: Comparison of experimental and theoretical static magnetic dipole (/x) in /in and electric quadrupole (Q) in e.fm 2 
moments of the ground states of the p-shell nuclei of interest. The quadrupole moments have been calculated either with 
harmonic oscillator (HO) or Woods-Saxon (WS) wave functions using method II with a volume term (see text), according to 
1151 1. The experimental data are taken from Ref. [4t|, if not specified. 
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TABLE II: Comparison of experimental and theoretical B(GT)-values in p-shell mirror nuclei. The results shown are obtained 
with the isospin-nonconserving (INC) interaction and either harmonic oscillator (HO) or Woods-Saxon (WS) wave functions 
determined with method I or II (see text) by adjusting the volume term, without including the dependence on the parent state. 
The results obtained by including a surface term are shown into parenthesis. The B(GT)-values given in column (INC+WS)* 
incorporate the dependence on the parent states as follows from Eq. JSJ. The experimental values [TBI are obtained following 
the procedure of Refs. |3Ct l31f on the basis of Q- values from Ref. |32|, halflives and branching ratios from Ref. |31f . if not 
specified. 
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TABLE III: Comparison of experimental and theoretical B(GT)-values in sd-shell mirror nuclei. We present results obtained 
with the isospin conserving (IC) and isospin non-conserving (INC) interactions and either harmonic oscillator (HO) or Woods- 
Saxon (WS) wave functions determined from method II by adjusting a volume term. The results including the dependence on 
the parent state are shown in column (INC+WS)* only for the nuclei considered in the discussion of the asymmetries. The 
experimental values fl5| are obtained following the procedure of Refs. |3(ll3l| on the basis of Q- values from Ref. |32l| . halflives 
and branching ratios as specified. 
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0.0052 


0.0463 
0.0830 
0.0643 
0.0047 






25 Si(/3+) 25 Al 


5 + 
2 


3 
2 


5 + 

2 1 

3 + 
2 1 
7 + 
2 1 
5 + 
22 


la 
i 

2 
1 
2 
1 

2 


;.s.) 


0.0463 
0.0765 
0.0589 
0.0049 


0.0462 
0.0798 
0.0595 
0.0053 


0.0458 
0.0738 
0.0578 
0.0050 






28 Al(/T) 28 Si 


3+ 


1 









0.1458 


0.1444 


0.1435 




0.0825(1) [37] 


28 P(/3+) 28 Si 


3+ 


1 


2+- 

Z l ! 







0.1449 


0.1303 


0.1287 




0.0870(9) [2JJ 


31 Al(/3-) 31 Si 


5 + 
2 


5 
2 


3 + 
2 1 


1 (g-s.) 


0.2120 


0.2090 


0.2124 




0.1049(122) [3J 


31 Ar(/3+) 31 Cl 


5 + 
2 


5 
2 


3 + 
2 1 


1 (g-s.) 


0.1947 


0.1958 


0.1827 




0.0788(53) [38] 


35 S(/T) 35 C1 


3 + 
2 


3 
2 


3 + 
2 1 


la 


v s.) 


0.0861 


0.0870 


0.0865 




0.0588(2) [37] 


35 K(/3+) 35 Ar 


3 + 
2 


3 
2 


3 + 
2 1 


la 


;.s.) 


0.0846 


0.0843 


0.0824 




0.0496(84) [37] 


35 P(/T) 35 S 


1 + 
2 


5 
2 


1 + 

2 1 


3 
2 




0.9052 


0.8572 


0.8836 


0.8861 


0.4591(70) [22] 


35 Ca(/3+) 35 K 


1 + 
2 


5 
2 


1 + 

2 1 


3 
2 




0.8572 


0.7435 


0.7257 


0.7458 


0.3891(154) [22] 
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TABLE IV: Comparison of experimental and theoretical asymmetries 8 (%) of B(GT)-values for transitions in selected p— 
and sd-shell mirror nuclei. Entries denoted A = 20, A = 20(a) and A = 35 correspond to O-F-Mg-Na, F-Ne-Na and P- 
S-Ca-K systems, respectively (see Table II for the details). The different columns give the asymmetries obtained with the 
isospin conserving (IC) or isospin non-conserving (INC) interactions and harmonic oscillator (HO) or Woods-Saxon (WS) wave 
functions obtained with method II (see text), adjusting a volume term. The column (INC+WS)* shows the results obtained 
with the inclusion of the dependence on t he p arent state. The experimental asymmetries are calculated on the basis of the 
data from Tables II - III (see also are Ref. |l5|p. 

A Jf;Ti JJ;T f INC+HO IC+WS INC+WS (INC+WS)* EXP. 



8 


2+ 


1 




) 




3.10 


13.21 


16.65 


18.82 


8.4+1.8 


9 


3 - 
2 


. 3 

' 2 


3- 
21 


1 

2 


(g-s.) 


-1.26 


6.4 


4.95 


2.72 


16 + 8 


12 


1+ 


1 


of; 


0( 


g-s.) 


0.70 


9.4 


10.09 


21.85 


12.6 + 0.8 








2x5 


) 




-0.65 


8.45 


7.92 


11.33 


9.5 + 1.9 


13 


3 - 
2 


. 3 
' 2 


1 - 

2 1 


1 

2 


(g- s -) 


2.10 


7.26 


9.50 


8.65 


11.3 + 3.2 


17 


1 - 
2 


. 3 

' 2 


3- 
2 1 


1 
2 




53.26 


5.29 


58.90 


64.80 


44 + 7 


20 


0+ 


2 


if; 


1 




0.04 


1.41 


1.42 


-1.53 


18 + 8 


20(a) 


2+ 


1 









-1.45 


-13.64 


-14.99 


1.11 


1.8 + 0.7 


35 


i + 

2 


■ ' r > 

' 2 


1 + 

2 1 


3 
2 




5.59 


15.28 


21.75 


18.81 


18 + 5 



TABLE V: Comparison of experimental and theoretical asymmetries S (%) of B(GT)-values for transitions in sd-shell mirror 
nuclei. The entry denoted as A = 35 correspond to the S-Cl-K-Ar system (see Table II for the details). The different 
columns give the asymmetries obtained with the isospin conserving (IC) and isospin non-conserving (INC) interactions with 
the harmonic oscillator (HO) or Woods-Saxon (WS) wave functions obtained with method II (see text), adjusting a volume 
term. The experimental asymmetries are calculated on the basis of the data from Tables II - III (see also are Ref. 15]). 



A Jf;Tj JJ;Tf INC+HO IC+WS INC+WS EXP. 



21 


5 4 

2 


. 3 

> 2 


3 + 
2 1 
5 + 
2 1 


1 

2 
1 
2 


(g-s-) 


-3.96 
-1.71 


3.49 
2.35 


-0.56 
0.60 








7 + 
2 1 


1 

2 




2.78 


1.44 


4.23 


25 


5 4 

2 


. 3 
■ 2 


5 4 
2 1 
M + 
2 1 


1 

2 
1 

2 


(g-s-) 


-1.09 
9.21 


2.24 
2.83 


1.11 
12.39 








7 + 
2 1 


1 

2 




7.76 


3.83 


11.23 








5 4 
22 


1 

2 




-3.17 


-2.52 


-5.58 



28 3+;l 2+;0 0.61 10.85 11.54 -5 + 1 

31 f + ;f § + ;§(g.s.) 8.88 6.72 16.27 33 + 18 
35| + ;| |+;i(g.s.) 1.82 3.13 4.96 18.4 + 20.0 
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FIG. 1: Experimental asymmetry J1J deduced from the data on /3± mirror transitions characterized by |AT| = 1, |AJ| = 0, 1, 
7Tj = 717 (calculated on the basis of the experimental B(GT)-values given in Tables II and III; see the captions of the tables for 
references). The results for the second transitions in A = 12, A — 20 O-F-Mg-Na and A = 35 S-Cl-K-Ar are shifted by 0.5 unit 
from the A- value. 




FIG. 2: Experimental versus theoretical (Table I): (a) absolute values of magnetic dipole moments of the ground states of 
p-shell nuclei obtained with the free-nucleon (/-factors (data are taken from Ref. |46l l48t l5Cjp: (b) absolute values of electric 
quadrupole moments of the ground states of p-shell nuclei obtained with the WS wave functions and e(n) = 1.35, e(y) = 0.35 
(data are taken from Ref. |4o. I47ll49ll5l|'l. Solid lines correspond to the ratios equal to 1 and are given to guide the eye. 
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FIG. 3: Ratios of experimental to theoretical B(GT) values present in Tables II and III (see the captions of the tables for 
references on the experimental data). Theoretical values are those obtained with the dependence on the parent state (column 
(INC+WS)* in Tables II and III) and without dependence for A = 28, 31 and A = 35 S-CL-K-Ar (column (INC+WS) in Table 
III) . Solid and dotted lines correspond to the rations 1 and (0.77) 2 , respectively, and are given to guide the eye. 




FIG. 4: Experimental versus theoretical asymmetries S (left) and ratios of experimental to theoretical asymmetries (right) of 
selected GT decay rates for p- and sd-shell nuclei. For a given mass, the theoretical points show the effect due to the inclusion 
of the dependence on the parent state, except for A = 28,31 and A = 35 S-Cl-K-Ar system. The results for the second 
transitions in A — 12 and A = 35 S-Cl-K-Ar are shifted by 0.5 unit from the A-value. The large in absolute value ratio of 
the asymmetries for A = 20 case O-F-Mg-Na is excluded. The experimental asymmetries are calculated on the basis of the 
experimental B(GT)-values from Tables II and III (see the captions of the tables for references). 



